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METHODS AND APPARATUS FOR ESTIMATING PHYSICAL PARAMETERS OF 
RESERVOIRS USING PRESSURE TRANSIENT FRACTURE INJECTION/FALLOFF 

TEST ANALYSIS 

Field of the Invention 



[0001] The present invention pertains generally to the field of oil and gas subsurface 

earth formation evaluation techniques and more particularly to a method and an apparatus for 
evaluating physical parameters of a reservoir using pressure transient fracture injection/falloff 
test analysis. More specifically, the invention relates to improved methods and apparatus using 
graphs of transformed pressure and time to estimate permeability and fracture-face resistance of 
a reservoir. 
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Background of the Invention 

[0002] The oil and gas products that are contained, for example, in sandstone earth 

formations, occupy pore spaces in the rock. The pore spaces are interconnected and have a 
certain permeability, which is a measure of the ability of the rock to transmit fluid flow. When 
some damage has been done to the formation material immediately surrounding the bore hole 
during the drilling process or if permeability is low, a hydraulic fracturing operation can be 
performed to increase the production from the well. Hydraulic fracturing is a process by which 
a fluid under high pressure is injected into the formation to split the rock and create fractures 
that penetrate deeply into the formation. These fractures create flow channels to improve the 
near term productivity of the well. 

[0003] Evaluating physical parameters of a reservoir play a key part in the appraisal of 

the quality of the reservoir. However, the delays linked with these types of measurements are 
often very long and thus incompatible with the reactivity required for the success of such 
appraisal developments. 

[0004] One of the reasons is the complexity of a multilayer environment, it increases as 

the number of layers with different properties increases. Layers with different pore pressure, 
fracture pressure, and permeability can coexist in the same group of layers. The biggest 
detriment for investigating layer properties is a lack of cost-effective diagnostics for 
determining layer permeability, and fracture-face resistance of reservoir. 
[0005] Numerous analyses have been carried out to evaluate physical parameters of a 

reservoir. More particularly, before-closure pressure-transient analysis has been commonly 
used to estimate permeability and fracture-face resistance from the pressure decline following a 
fracture-injection/falloff test in the reservoir. 

[0006] Before-closure pressure-transient analysis is described by Mayerhofer and 

Economides in a paper SPE 26039 "Permeability Estimation From Fracture Calibration 
Treatments," presented at the 1993 Western Regional Meeting, Anchorage, Alaska, 26- 
28 May 1993; also by Mayerhofer, Ehlig-Economides, and Economides in a journal JPT 
(March 1995) on page 229 "Pressure-Transient Analysis of Fracture-Calibration Tests"; and by 
Ehlig-Economides, Fan, and Economides in a paper SPE 28690 "Interpretation Model for 
Fracture Calibration Tests in Naturally Fractured Reservoirs" presented at the 1994 SPE 
International Petroleum Conference and Exhibition of Mexico, 10-13 October 1994 . The 
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analysis was formulated in part using the early-time infinite-conductivity fracture solution of 
the partial differential equation that Gringarten, Ramey, and Raghavan suggested in a journal 
SPEJ (August 1974) on page 347 "Unsteady-State Pressure Distributions Created by a Well 
With a Single Infinite-Conductivity Vertical Fracture" which assumed the use of a slightly 
compressible reservoir fluid. However, diagnostic fracture-injection/falloff tests are commonly 
implemented in reservoirs containing highly compressible fluids, for example, in natural gas 
reservoirs. When the compressibility of the reservoir fluid deviates from the assumption of a 
slightly compressible fluid, the analysis methods as used in the prior art can lead to erroneous 
permeability and fracture-face resistance estimates. 

[0007] The errors in the estimates of the permeability and fracture-face resistance are 

significant and can be detected in the plotting of the experimental data obtained with a slightly 
compressible reservoir fluid. As a matter of fact, these errors are the consequences of the 
inaccuracy of the approximations as used in the prior art. These approximations used in 
connection with the actual theory developed with the pressure-transient leakoff analysis are 
based on the assumption that the reservoir fluid properties are not functions of pressure, which 
could not be the case when the reservoir fluid is a gas. The approximations as assumed in the 
prior art are as follows: 

1) Before-Closure Pressure-Transient Leakoff' Analysis Assuming a Slightly-Compressible 
Reservoir Fluid 

[0008] The pressure decline following a fracture-injection/falloff test can be divided 

into two distinct regions: before-fracture closure and after-fracture closure. Before-closure 
pressure-transient analysis is used to determine permeability from the before-fracture closure 
decline data. Mayerhofer and Economides in paper SPE 26039 divide the before-closure 
pressure difference between a point in an open, infinite-conductivity fracture and a point in the 
undisturbed reservoir into four components written as: 

Ap(t) = Ap res (0 + Ap cake (0 + Ap piz (t) + Ap fiz (t) ( 1 ) 

[0009] The pressure difference in the polymer invaded zom,&p piz (t), the filtrate 

invaded zone, Apf lz (t), and across the filtercake, Ap cake (t\ can be grouped into a fracture-face 
pressure difference term, Apf ace (t). Consequently, the pressure gradient consists of reservoir 
and fracture-face resistance components, and is written as: 

AP (0 = APres (0 + AP/oce (0 (2) 



2 



2) Fracture-Face Pressure Difference 

[0010] In the same way, in paper SPE 26039 Mayerhofer and Economides determine 

the fracture-face resistance pressure difference by using the concept of a fracture-face skin 
proposed by Cinco-Ley and Samaniego in paper SPE 10179 "Transient Pressure Analysis: 
Finite Conductivity Fracture Case Versus Damage Fracture Case" presented at the 1981 SPE 
Annual Technical Conference and Exhibition, San Antonio, Texas, 5-7 October 1981. Cinco- 
Ley and Samaniego defined fracture-face skin as: 



nb fs 



k fs 



•0) 



2L f 

where bf s is the damaged zone width, Lf is the fracture half-length, k is the reservoir 
permeability, andfy 5 is the damaged-zone permeability. Mayerhofer and Economides account 
for variable fracture-face skin by defining resistance, in paper SPE 26039, as: 

bf s (t) 

¥0=f^. (4) 

K fs 

and dimensionless resistance in journal JPT of (March 1995) by: 



(5) 



where R$ is the reference filtercake resistance at the end of the injection and ^ is the time at the 
end of the injection. 

[001 1] With Eqs. 4 and 5, fracture-face skin is written as: 

or as: 

[0012] Fracture- face skin is equivalent to a dimensionless pressure difference across the 

fracture face; thus, it can be written as: 
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"to-w^r 1 (8) 

where h p is the permeable reservoir thickness, q L ^ is the total injection (leakoff) rate into both 

wings of the hydraulic fracture, B is the formation volume factor of the filtrate, and /i is the 
filtrate viscosity. With Eq. 8, the fracture-face pressure difference is written as: 

^.MuwjaLffciE ,9) 

h p Lf I \t ne 

[0013] With a fracture symmetric about the wellbore, the total injection (leakoff) rate 

can be written as: 

U f B = 2q e (10) 

where ^ is the leakoff rate in one wing of the fracture. The fracture-face pressure difference is 
written as: 

t# face =U\2(n)-£%-q t JZ (11) 

h p L f Vne 

[0014] Define: 

Ro=mRq (12) 

where Rq is the fracture- face resistance, then the fracture-face pressure difference is written as: 

Ap face =l4L2^)-^-q e [±. (13) 

[0015] Assuming the fracture-face skin is a steady-state skin, the pressure difference at 

the fracture face at any time since the injection began is written as: 

(*Pface)n = 141 .2(^)-^-(^) n jX (14) 

h p L f V 'ne 

where the subscript n denotes a time /„ . 
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[0016] According to Nolte, K.G. in a journal SPEFE (December 1986): "A General 

Analysis of Fracturing Pressure Decline With Application to Three Models," on page 571, the 
leakoff rate from one wing of a hydraulic fracture during a shut-in period is written as: 

(qg) .=Jj±j\±\f!m} jj±_-\±ipj-i-pj) (15) 

Ht J [5.615 jSf L</(A0j - L5.615J ^ (tj-tj-i) V ' 

whereby is the fracture area, Sy is the fracture stiffness and the subscript y is a time index. Sf 

can be determined using Table 1 which summarizes what Valko and Economides determine in 
Chap. 2, pages 19-51: "Linear Elasticity, Fracture Shapes, and Induced Stresses," Hydraulic 
Fracture Mechanics, John Wiley & Sons, New York City (1997). The fracture stiffness Sf for 

2D fracture models can be calculated by using either one of the three formulas as shown in 
Table 1, the radial equation, the Perkins-Kern-Nordgren equation, or the Geertsma-deKlerk 
equation. 

[0017] Define: 
_ (Pj-i-Pj) 

J -^J^o (16) 

then the leakoff rate from one wing can be written as: 



^-jztst;'' < 17 > 



24 Af 

[0018] At any time during the shut-in period, t n >t ne , the fracture-face pressure 

difference is written as: 



V#face)n- hpLf JJ d n^~ (18) 



[0019] The ratio of permeable fracture area to total fracture area is defined by: 

A, 



.(19) 



where for a rectangular-shaped fracture,^ =h p Lf y md the fracture-face pressure difference at 
any time during the shut-in period, t n > t ne , is written as: 
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J 5.615 r p Sf \t ne 

[0020] Eq. 20 is also applicable to radial, elliptical, or other idealized fracture geometry 

by defining fracture- face skin in terms of equivalent fracture half-length, L e , and noting that any 
fracture area can be expressed in terms of an "equivalent" rectangular fracture area. 

3) Reservoir Pressure Difference 

[0021] As in previously mentioned article of the journal SPEJ (August 1974) on page 

347: "Unsteady-State Pressure Distributions Created by a Well With a Single Infinite- 
Conductivity Vertical Fracture", the pressure drop in the reservoir is modeled by Gringarten, 
Ramey, and Raghavan for a slightly-compressible fluid, and is written in dimensionless form 
as: 

PL f D = ] nt i f D . (21) 

where 



L f U ~\A\2q Lf B M 



and 



*LfD =0.0002637 ~~~y~ (23) 

<f>Hc t LJ 

[0022] In Eq. 23, ^ is the porosity andc, is the total compressibility. Equating Eqs. 21 

and 22 and combining with Eq. 10 results in: 

BA Pres =141.2(2)^- %> pZ^ (24) 

[0023] By expanding the dimensionless time term, the reservoir pressure difference can 

be written as: 

&p res = 141.2(2X0.02878) 1 Tf jj-ft^ (25) 

hpLfyjk \<pc t 

[0024] The pressure difference at any time t n is written using superposition as: 
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(Ap«,)„ = 141.2(2)(0.02878) / £f [(^). -{q t ) H l[i~T~ x 



.(26) 



[0025] In a simplification of the more general method, Mayerhofer and Economides in 

paper SPE 26039, and Valko and Economides in a journal SPEPF (May 1999) on page 117: 
"Fluid-Leakoff Delineation in High-Permeability Fracturing M , assume that during the injection, 
the first + 1 leakoff rates are constant, where ne is the index corresponding to the time at the 
end of the injection and the beginning of the pressure falloff, the leakoff rates can be written as: 

(q e )j = Constant 1 < j < ne + 1 , and (qi)o = 0 (27) 

[0026] With Eq. 27, the reservoir pressure difference at any time t n is written as: 



(Ap res ) n =141.2(2)(0.02878) 



or written as: 



h p L f 4k\<l>c t 



+ Z t(Vt)j-tot)j-\^'n-tj-\ 
j=ne+3 



(28) 



(Ap res ) n = 141.2(2)(0.02878) 



h p Lf4k \<pc t 



n *ne+\ 



+ Z Mj-teOj-l^'n-tj-l 

/ i 

*ne+\ 



+(%Wi>/£fi-J^f 



.(28) 



[0027] With Eq. 17 substituted for leakoff rate and Eq. 19 for the ratio of permeable to 

total fracture area, the reservoir pressure difference at any time^ is written as: 



_ 141. 2(2X0.02878)(24) 1_ [J 

5.615 r p Sfy]k\<pc t 



+ Z Wj-dj-l^n-tj-l 
j=ne+2> 



\ V l n J 

4) Specialized Cartesian Graph for Determining Permeability and Fracture-Face Resistance 



.(29) 
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[0028] Eq. 2 defines the total pressure difference between a point in the fracture and a 

point in the undisturbed reservoir as the sum of the reservoir and fracture-face pressure 
differences, which is written as: 



(4P)„ = 



141.2(2)(0 t Q2878)(24) 1_ 

5.615 r p Sf4k^<pc t 



n *ne+\ 



+ Z ldj-<tj-\]<J*n-tj-\ 



141.2(^)24 Rq 



5.615 r n S 



P f V l ne 



(30) 
[0029] 



Algebraic manipulation allows Eq. 30 to be written as: 



(Ap) n _141.2(2)(0.02878)(24) 1 



dnVn^e 5615 r p S f Jk\<t>c t 



dne+2 ( *n ~ l ne+l 1 
d n I Wne ) 



1/2 



j=ne+2> n 



'n-tj-l 
V *r£ne ) 



1/2 



7 ne+\ 



l-.i 



t ne+\ 



141.2(^)24 flp 1 
5.615 r p Sft m 



(31) 
[0030] 



In view of Eq. 16, the term</„ e+1 can be written in an alternative form as: 



_ 5.615 Sf 24 Af _ 5.615 Sf 



2 ne+\ ~ ' 



24 A f 5.615 Sf 



24 A 



(32) 



but recognizing that q ne = q ne+x andV Lne =(ge)„ e t ne allows Eq. 32 to be written as: 



.(33) 



24 < wc /!/ 

where v Lne is the leakoff volume at the end of the injection. Define lost width due to leakoff at 
the end of the injection as: 



and Eq. 33 can be written as: 



•(34) 
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, 5.615 1 
^\=—S f y> L - 



(35) 



ne 



[0031] 



Define: 



(36) 

(37) 
.(38) 



X„ = 



d ne+2 



V/2 



y=ne+3 



/ -/• .V /2 



\ *n*ne ) 



c 2 



d n'ne v 



1-. 1- 



.(39) 



141.2(2)(0.02878)(24) 1 

m M ~ — 



and 



[0032] 



5.615 



141.2(^)24 /fr 1 



Sfyfk 1 



r p *f 



.(40) 



5.615 r p Sf t ne 

Combining Eq. 31 and Eqs. 36 through 41 results in: 



.(41) 



y n = ™M x n+ b M (42) 

[0033] Eq. 42 suggests a graph of y n versus x n using the observed fracture- 

injection/falloff before-closure data will result in a straight line with the slope a function of 
permeability and the intercept a function of fracture-face resistance. Eqs. 41 and 42 are used to 
determine permeability and fracture-face resistance from the slope and intercept of a straight- 
line through the observed data. 



5) Before-Closure Pressure-Transient Leakoff Analysis in a Dual-Porosity Reservoir System 
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[0034] In the present application, dual porosity refers to a mathematical model of a 

naturally fractured reservoir system. In paper SPE 28690, Ehlig-Economides, Fan, and 
Economides formulated the Mayerhofer and Economides model for dual-porosity reservoirs 
using Cinco-Ley and Meng's dimensionless pressure. In a paper SPE 18172: "Pressure 
Transient Analysis of Wells With Finite Conductivity Vertical Fractures in Dual Porosity 
Reservoirs," presented at the 1988 SPE Annual Technical Conference and Exhibition, Houston, 
Texas, 2-5 October 1988, Cinco-Ley and Meng determine dimensionless pressure with an 
early-time approximation for flow of a slightly compressible fluid from an infinite-conductivity 
fracture as: 



\ nt LrD 

>h°=\-l (43) 

where for dual-porosity reservoirs, 

k jb h p *Pres .... 

t L £,=0.0002637-^- (45) 

4>Hc t Lj 

and<yis the natural fracture storativity ratio as defined by Warren, J.E. and Root, P.J. in a 
journal SPE J (September 1963) on page 245: "The Behavior of Naturally Fractured 
Reservoirs". 

[0035] Writing Eq. 43 as 

°>PL f D = ^L f D (46) 

and repeating the derivation for the reservoir pressure difference results in changing the final 
slope definition, Eq. 40, to: 

_ 141.2(2)(0.02878)(24) 1 ... 
5.615 rpSf^Jakp 

[0036] In a dual-porosity reservoir or in a naturally fractured reservoir system, before- 

closure pressure-transient leakoff analysis using the specialized Cartesian graph results in an 
estimate oiakp,. Methods as used in the prior art allow the product to be evaluated without an 



10 



acceptable accuracy, and estimating fracture storativity co or bulk-fracture permeability kp, 
requires additional testing which would involve additional inaccuracy. Therefore, since the 
permeability and fracture-face resistance evaluations cannot be directly obtained and since the 
additional testing increase the error of these evaluations, it is necessary to determine the 
product cokjb with more accuracy. 

[0037] Henceforth, there is a need to find another approach that mitigates nonideal 

leakoff behavior attributed to pressure-dependent fluid properties with more accuracy. For 
example, in low pressure gas reservoirs, that is, in many gas reservoirs with a pore pressure less 
than about 3000 psi, reservoir fluid properties are strong functions of pressure. When fluid 
properties are strong functions of pressure, assuming constant properties for use in pressure and 
time formulations will cause significant error in permeability and fracture-face resistance 
determinations. 

[0038] These approximations as used in the prior art are therefore unsatisfactory. Thus, 

there is a desire not only for estimating accurate permeability and fracture-face resistance of a 
reservoir to appraise its quality but also for avoiding the delays linked with this type of 
measurements which are often very long and incompatible with the reactivity required for the 
success of such appraisal developments. New, faster and accurate evaluation means are 
therefore sought as a decision-making support. 
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Summary of the Invention 

[0039] The present invention pertains to a method and an apparatus for evaluating 

physical parameters of a reservoir using pressure transient fracture injection/falloff test 
analysis. 

[0040] The before-closure pressure-transient leakoff analysis for a fracture- 

injection/falloff test is used to mitigate the detrimental effects of pressure-dependent fluid 
properties on the evaluation of the permeability and fracture-face resistance of a reservoir. A 
fracture-injection/falloff test consists of an injection of liquid, gas, or a combination (foam, 
emulsion, etc.) containing desirable additives for compatibility with the formation at an 
injection pressure exceeding the formation fracture pressure followed by a shut-in period. The 
pressure falloff during the shut-in period is measured and analyzed to determine permeability 
and fracture-face resistance by preparing a specialized Cartesian graph from the shut-in data 
using adjusted pseudovariables such as adjusted pseudopressure data and adjusted pseudotime 
data. This analysis allows the data on the graph to fall along a straight line with either constant 
or pressure-dependent fluid properties. The slope and the intercept of the straight line are 
respectively indicative of the permeability and fracture-face resistance evaluations. 
[0041] Pseudo variable formulations for before-closure pressure-transient fracture- 

injection/falloff test analysis minimize error associated with pressure-dependent fluid 
properties by removing the "nonlinearity". The use of adjusted pseudovariables according to the 
present invention allows analysis to be carried out when a compressible or slightly 
compressible fluid is injected into a reservoir containing a compressible fluid. Therefore, the 
permeability and the fracture-face resistance of the reservoir can be estimated with more 
accuracy by the pressure transient fracture injection/falloff test. 

[0042] Although the primary benefit occurs when the reservoir fluid is highly 

compressible, the technique is also valid for all reservoir fluids that are either compressible or 
slightly compressible. 

[0043] In accordance with a first aspect of the present invention, a method of estimating 

physical parameters of porous rocks of a subterranean formation containing a compressible 
reservoir fluid comprising the steps of injecting an injection fluid into the subterranean 
formation at an injection pressure exceeding the subterranean formation fracture pressure, 
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shutting in the subterranean formation, gathering pressure measurement data over time from the 
subterranean formation during shut-in, transforming the pressure measurement data into 
corresponding adjusted pseudopressure data to minimize error associated with pressure- 
dependent reservoir fluid properties, and determining the physical parameters of the 
subterranean formation from the adjusted pseudopressure data. 

[0044] In an embodiment, the adjusted pseudopressure data is defined by the equation: 

*(Pw)n 



P 



Pdp 



[0045] Furthermore, the determination of the physical parameters is obtained by a plot 

of the adjusted pseudopressure data over time showing a straight line characterized by a slope 
m M and an intercept b M , wherein % is a function of permeability k and is a function of 
fracture-face resistance Ro wherein: 

k = 



(141.2)(2)(0.02878)(24) 1 



5.615 r p Sfm M 
D 5615 

*° = 141.2*(24) r P S f^ bM - 

[0046] In accordance with a second aspect of the present invention, a method of 

estimating physical parameters of porous rocks of a subterranean formation containing a 
compressible reservoir fluid comprising the steps of injecting an injection fluid into the 
subterranean formation at an injection pressure exceeding the subterranean formation fracture 
pressure, shutting in the subterranean formation, gathering pressure measurement data over 
time from the subterranean formation during shut-in, transforming the pressure measurement 
data into corresponding adjusted pseudopressure data and time into adjusted pseudotime data to 
minimize error associated with pressure-dependent reservoir fluid properties, and determining 
the physical parameters of the subterranean formation from the adjusted pseudopressure data. 
[0047] In an embodiment, the adjusted pseudopressure data and the adjusted 

pseudotime are defined by the equations: 



^ dAt , 
, and 

0 (Mg c t)w 
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[0048] w ,=Mri* 

P Jo M g c, 



[0049] Furthermore, the determination of the physical parameters is obtained by a plot 

of the adjusted pseudopressure data over adjusted pseudotime data showing a straight line 
characterized by a slope m M and an intercept b M , wherein m M is a function of permeability k and 
bM is a function of fracture-face resistance R 0 wherein: 

(141.2)(2)(0.02878)(24) 1 



5.615 r p Sftn M 
5.615 

^ S 141.2»(24)* A ' 

[0050] Also in one embodiment, the reservoir fluid is compressible or slightly 

compressible. 

[0051] And in another embodiment, the injection fluid is compressible or slightly 

compressible. 

[0052] In accordance with a third aspect of the present invention, a system for 

estimating physical parameters of porous rocks of a subterranean formation containing a 
compressible reservoir fluid comprising a pump for injecting an injection fluid into the 
subterranean formation at an injection pressure exceeding the subterranean formation fracture 
pressure, means for gathering pressure measurement data from the subterranean formation 
during a shut-in period, means for transforming the pressure measurement data into adjusted 
pseudopressure data to minimize error associated with pressure-dependent reservoir fluid 
properties and means for determining the physical parameters of the subterranean formation 
from the adjusted pseudopressure data. 

[0053] In an embodiment, the determining means comprises graphics means for 

plotting a graph of the adjusted pseudopressure data over time, the graph representing a straight 
line with a slope m M and an intercept b M wherein m^is a function of permeability k and b M is a 
function of fracture-face resistance R 0 . 

[0054] In accordance with a fourth aspect of the present invention, a system for 

estimating physical parameters of porous rocks of a subterranean formation containing a 
compressible reservoir fluid comprising a pump for injecting an injection fluid into the 
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subterranean formation at an injection pressure exceeding the subterranean formation fracture 
pressure, means for gathering pressure measurement data from the subterranean formation 
during a shut-in period, means for transforming the pressure measurement data into adjusted 
pseudopressure data and time into adjusted pseudotime to minimize error associated with 
pressure-dependent reservoir fluid properties and means for determining the physical 
parameters of the subterranean formation from the adjusted pseudopressure data. 
[0055] In an embodiment, the determining means comprises graphics means for 

plotting a graph of the adjusted pseudopressure data over adjusted pseudotime data, the graph 
representing a straight line with a slope niM and an intercept wherein is a function of 
permeability k and bu is a function of fracture-face resistance Ro. 

[0056] Also in another embodiment, the reservoir fluid is compressible or slightly 

compressible. 

[0057] And in another embodiment, the injection fluid is compressible or slightly 

compressible. 

[0058] Other aspects and features of the invention will become apparent from 

consideration of the following detailed description taken in conjunction with the accompanying 
drawings. 
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Brief Description of the Drawings 

[0059] A more complete understanding of the present disclosure and advantages thereof 

may be acquired by referring to the following description taken in conjunction with the 
accompanying drawings wherein: 

[0060] Figure 1 shows a Tablel representing three formulas used for the calculation of 

fracture stiffness for 2D fracture models. 

[0061] Figure 2 shows a Table2A which lists equations and definitions for before- 

closure pressure-transient fracture injection/falloff test analysis. 

[0062] Figure 3 shows a Table2B which lists additional equations and definitions for 

before-closure pressure-transient fracture injection/falloff test analysis. 

[0063] Figure 4 shows a plotting of three specialized Cartesian graphs of the basic 

linear equations y n versus x n according to a first series of experiments. 

[0064] Figure 5 shows a plotting of three specialized Cartesian graphs of the basic 

linear equations y n versus x n according to a second series of experiments. 

i- 

[0065] Figures 6A, 6B and 6C are a general flow chart representing a method of 

iterating the measurements and plotting the Cartesian graphs thereof. 

[0066] Figure 7 shows schematically an apparatus located in a wellbore useful in 

performing the methods of the present invention. 

[0067] The present invention may be susceptible to various modifications and 

alternative forms. Specific embodiments of the present invention are shown by way of 
example in the drawings and are described herein in detail. It should be understood, however, 
that the description set forth herein of specific embodiments is not intended to limit the present 
invention to the particular forms disclosed. Rather, all modifications, alternatives and 
equivalents falling within the spirit and scope of the invention as defined by the appended 
claims are intended to be covered. 
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Description of The Embodiments of The Invention 

[0068] The methods as shown in the prior art for analyzing the before-closure pressure 

decline following a fracture-injection/falloff test do not consider a compressible reservoir fluid 
with either a slightly compressible or compressible injection fluid. Accounting for 
compressible fluids is accomplished by using pseudovariables, or for convenience, adjusted 
pseudovariables in the derivation. 

[0069] Pseudovariables have been demonstrated in other well testing applications as 

removing the "nonlinearity" associated with pressure-dependent fluid properties, and using 
pseudovariable formulations for before-closure pressure-transient fracture-injection/falloff test 
analysis will minimize error associated with pressure-dependent fluid properties. Definitions of 
pseudovariables and adjusted pseudovariables can respectively be found in a paper SPE 8279 
by Agarwal, R.G.: "Real Gas Pseudo-time — A New Function for Pressure Buildup Analysis of 
MHF Gas Wells" presented at the 1979 SPE Annual Fall Technical Conference and Exhibition, 
Las Vegas, Nevada, 23-26 September 1979, and in a journal PEFE (December 1987) on page 
629 by Meunier, D.F., Kabir, C.S., and Wittman, M.J.: "Gas Well Test Analysis: Use of 
Normalized Pseudovariables". 

[0070] As a matter of fact, since Gas viscosity, deviation factor (z), and compressibility 

are functions of pressure; thus the governing partial differential equation is nonlinear. 
Therefore, pseudopressure and pseudotime are required to linearize the partial differential 
equation corresponding to the solution that Gringarten, Ramey, and Raghavan suggested in 
previously mentioned journal SPEJ (August 1974). Pseudopressure "corrects" for gas viscosity 
and real-gas deviation factor, and pseudotime "corrects" for gas viscosity and gas 
compressibility. Some authors find the use of pseudotime unnecessary as gas compressibility is 
nearly constant in most applications; however, both pseudopressure and pseudotime must be 
used to rigorously transform the governing partial differential equation to a linear partial 
differential equation. 

[0071] Using both pseudopressure and pseudotime enables well design engineers to 

obtain the best "correct" answer. However acceptable answers may be obtained using only 
pseudopressure. Two series of experiment will be shown later in Figures 4 and 5 which 



17 



illustrate three graphs resulting in the evaluation of permeability and fracture-face resistance 
when pressure and time; pseudopressure and time; and finally pseudopressure and pseudotime 
formulations represent the variables. 

1) Reservoir Adjusted Pseudopressure Variables Difference 

[0072] For convenience, the new approach is illustrated with adjusted pseudovariables. 

The pressure drop in the reservoir modeled by Gringarten, Ramey, and Raghavan in SPEJ 
(August 1974) for a slightly-compressible fluid, is written in dimensionless form as: 

PL f D = ^ nt L f D • ( 48 )> th e same as Eq. 21. 
[0073] Writing Eq. 48 in terms of pseudopressure accounts for the variation of viscosity 

and gas deviation factor for the compressible fluid in the reservoir. Define adjusted 
pseudopressure variable as: 

p a -m p ^ (49) 

P Jo V g z 

where z is the gas deviation factor, JI is the viscosity evaluated at average reservoir pressure, z is 
the gas deviation factor at average reservoir pressure, and p is average reservoir pressure. The 
derivative of Eq. 49 is written as: 

^Pa_ = Ml_P_ = pB_^^Pa^ (50) 
dp p pi fi B Ap 

[0074] With Eq. 50, the definition of dimensionless pressure is written as: 

khpty _ kh p A Pa 

which when combined with Eq. 48 results in: 

PaLfD = yjxtLfD (52) 

[0075] The reservoir pressure difference in terms of adjusted pseudopressure variable 

can now be written as: 
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.(53) 



[0076] 
as: 



With Eq. 10, the reservoir adjusted pseudopressure variable difference is written 



(54) 



[0077] Dimensionless time is evaluated at average reservoir pressure, that is, 

dimensionless time is written as: 



*LfD 



0.00Q2637fo 
<t>p g c t I?f 



.(55) 



and the reservoir adjusted pseudopressure variable difference is written as: 



(Ap a ) res = 141.2(2)(0.02878)- 



1 



.(56) 



h p L f Jk\<t>c t B g 

[0078] The reservoir adjusted pseudopressure variable difference at any time t n is 

written using superposition as: 



[(*Pa)res]n = 141.2(2)(0.02878) 



.(57) 



[0079] The Valko and Economides assumption, in SPEPF (May 1999), that the 

first ne+1 leakoff rates are constant is modified such that the first ne+\ leakoff rates are constant 
at standard conditions. The assumption can now be expressed as: 



= Constant 1< j < ne + \, 



.(58) 



and implies that the pressure in the fracture during the injection is approximately constant. 
With Eq. 58, the reservoir adjusted pseudopressure variable difference at any time t n is written 
as: 
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K4P fl )«,]» = 14l.2(2)(0.02878) 



B, 



•g fg 



h p LfJky<t>c t 



(torts) 



n 

+ I 



B 



g J 



ne+2 



((«>*] 












I »g J 





f toth 

ft 

\ 8 Jne+\ 



(59) 



or: 



[(Ap a ) res ] w =141.2(2)(0.02878) 



B, 



8 \H 



h p Lf4k\<l>c t 



(H)g) 



ft 

^ 8 Jne+2 



j=ne+2> 



yfin *ne+\ 







(tot)g' 




{ h J 


j 


[ *g J 


H 



>g 



\ h j 



ne+\ 



1-. 1- 



,(60) 



[0080] The leakoff rate shown in Eq. 15 must be expressed in terms of adjusted 

pseudopressure variable, and is written as: 



[0081] 



24 ~\ Af (Mg B g)j (Pa)j-\ -(Pa)j 
5.615 ]Sf P g B g tj-t H 

Define: 



.(61) 



\ a atj - - ~ ' 

then Eq. 61 can be written as: 
[<%),]/- 



•(62) 



24 



5.615 



•(63) 



V °g 

[0082] With Eq. 63, the reservoir adjusted pseudopressure variable difference at any 

time?„ is written using superposition as: 
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K&Pa)res]n = 



141.2(2)(0.02878)(24) Af 1 \fg 
5.615 h p Lf Sfjk'i<l>c t 



{ d a ) ne+ 2 yl'n ~ l ne+\ 



j=ne+3 



t„\\-\\- 



.(64) 



or with Eq. 19, written as: 



141.2(2X0.02878)(24) 1_ \Mg_ 

5.615 r p Sf4ki</>c t 



(da ) ne +2 yfin 'ne+\ 

+ I [K)y-K) y _i]V'«-'y-i 

j=ne+3 



•(65) 



2) Fracture-Face Adjusted Pseudopressure Variable Difference 

[0083] The fracture-face adjusted pseudopressure variable difference is developed 

beginning from Eq. 8, which is written in terms of adjusted pseudopressure variable as: 



kh p(&Pa)face 
PaLrD = TTTZ; : H - = s f . 



.(66) 



or: 



(APoW=141 , w ^^rr 

h p L f 2 y\t m 

[0084] With Eq. 10 written for gas, the fracture-face adjusted pseudopressure variable 

difference is written as: 



.(67) 



(Ap fl )/ ace =141.2(^)^-f 



(68) 



and assuming a steady-state fracture-face skin, written as: 



K4P a )/fl«ln=141.20r)3£^ 



5, 



* J 



•(69) 
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for any timer,,. 



[0085] 



Define: 



*o»P*b (70) 

and the fracture-face adjusted pseudopressure variable difference is written as: 



B s Rq 
h p L f 



.(71) 



[0086] With Eq. 60 for the leakoff rate in terms of adjusted pseudopressure variable, the 

fracture-face adjusted pseudopressure variable difference is written as: 



K4p fl )/«*]#! 



_ 141.200(24) A f Rp 



(d a )n 



5.615 h p LfSf H "V'i 



.(72) 



or 



5.615 r p S f \\t m 

3) Specialized Cartesian Graph for Determining Permeability and Fracture-Face Resistance In 
Terms of Adjusted Pseudopressure Variable 

[0087] Eq. 2 defines the total pressure difference between a point in the fracture and a 

point in the undisturbed reservoir as the sum of the reservoir and fracture-face pressure 
differences, which is written in terms of adjusted pseudopressure variable as: 

&p a (t) = (Ap a ) res (t) + ( Ap a ) face (0 (74) 

[0088] Combining Eqs. 65, 73, and 74 results in the adjusted pseudopressure variable 

difference at any time t n , which is written as: 
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141.2(2)(0.02878)(24) 1 I Mg 
5.615 r p S f 4ki<fic t 



( d a ) ne +2 V'« ~ l ne+\ 

t [Wj-Mj_ x ]j^rx 



j=ne+7> 



141.2(^)(24) Rp 

5.615 r p s/ a)n 



(75) 
[0089] 



Algebraic manipulation of Eq. 75 results in: 



(Ap a )„ _141.2(2)(0.02878)(24) 1 Mg 



(d a ) n ^^ 5.615 



141.2(^)(24) Rq 1 



0 Sf4k\<l>c 



n 



1/2 



+ z 

j=ne+3 _ 



%-o-iV /2 



1-1- 



*ne+\ 



5.615 ^ 

[0090] The term (tf a ) w _+i can be written in an alternative form as: 



.(76) 



_ 5.615 S f B g 2 4 A f (B g ) ne + } 5 . 615 S f B g 

24 (B g ) ne+l 5.615 £ g 24 A f (B g ) ne+X * 

but recognizing thzt[(q e ) g / B] ne =[(qt) g / B] ne + { mdV Lne =[(qt) g ] ne t ne allows Eq. 77 to be written 
as: 



.(78) 



_5.615.-S-/ lg F_ ne 

where V Lne is the leakoff volume at the end of the injection. Define lost width due to leakoff at 
the end of the injection as: 



ne 



.(79) 



and Eq. 78 can be written as: 
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( d a)ne+\ = 

[0091] 



5.615 



8 



1 



24 ' "(B g ) ne t ne 
Define: 



(80) 



5.615 _ " g fi g 



(y a )n = 



24 J {B g ) ne ^<t>c, 



( V/2 



.(81) 

•(82) 
.(83) 



c al 



y'=ne+3 
c a2 



x .3/2 
\ a a'n l ne 



1- 



(84) 



and recall: 



141.2(2)(0.02878)(24) 1 

m M = — — 



5.615 



Sfjk ' 



(85) 



and 



141.2(^)24 Rq 



.(86) 



5.615 r p S f t ne 

[0092] Combining Eq. 76 and Eqs. 80 through 86 results in: 

(y a )n=m M (x a ) n +b M (87) 

[0093] Eq. 42 suggests a graph of(y a ) n versus (x a ) n using the observed fracture- 

injection/falloff before-closure data will result in a straight line with the slope a function of 
permeability and the intercept a function of fracture-face resistance, keeping in mind that the 
formulations of the slope does not change with the use of pseudovariables such that m M , m aM 
and m ap M are the same, but the values of the slope will change using the transformed pressure 



24 



measurement data. Eqs. 86 and 87 are used to determine permeability and fracture-face 
resistance from the slope and intercept of a straight-line through the observed data. 



4) Bef ore-Closure Pressure-Transient Leakoff Analysis in Dual-Porosity Reservoirs in Terms 
of Adjusted Pseudopressure Variable 

[0094] In dual-porosity reservoir systems, before-closure pressure-transient analysis in 

terms of adjusted pseudopressure variable changes by only one equation from the single- 
porosity case. Eq. 85 is modified and written as 

_ 141.2(2)(0.Q2878)(24) 1 (m 



5.615 



r, 



[0095] Consequently, the product of natural fracture storativity and bulk fracture 

permeability is determined from before-closure pressure-transient leakoff analysis in terms of 
adjusted pseudopressure variable in dual-porosity reservoir systems. 

[0096] A similar derivation can be used to derive the equations written in terms of 

adjusted pseudopressure and adjusted pseudotime variables. A similar derivation could also be 
used to demonstrate that other before-closure pressure transient analysis formulations can be 
expressed in terms of pseudovariables, but since most of the steps are the same, it would be 
redundant to repeat each derivation. 

[0097] Table2A in Figure 2 defines the parameters and variables used in the linear 

equations y n versus x n required for preparing the specialized Cartesian graphs in terms of 
pressure and time on a first column 212; adjusted pseudopressure variable and time on a second 
column 213; and adjusted pseudopressure and adjusted pseudotime variables on a third column 
214. 

[0098] For each of the three columns, pressure and time 212, adjusted pseudopressure 

variable and time 213, adjusted pseudopressure and adjusted pseudotime variables 214, the 
coefficients corresponding to the basic straight line equations are defined. These basic 
equations as shown in row 201, are respectively: y n =bj^ + m M x n , (y a )n =b M + m Af( x a)n 
and (y ap )n = b M + m M( x ap)n • 
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[0099] In the second row 202, the formulas of the before-closure pressure-transient 

analysis variable and adjusted variable with time and adjusted pseudotime variable 
yn> (ya)n> or (yap)n are respectively given as function of pressure p, pressure in reservoir p r and 

adjusted pseudopressure variable p Q and p ary and time at time step t n and at the end of an 
injection t ne . 

[00100] In the same way, in the third row 203, the formulas of the before-closure 
pressure-transient analysis variables and adjusted variables with time and adjusted pseudotime 
variable x n , (x a ) n , or (x ap ) n are respectively given as functions of coefficients (d a ) y (d ap \ (cy), 

(cai), (c ap i), (q), (c a2 ), {c ap 2) at time step t n9 at the end of an injection t ney or at the end of an 
adjusted pseudotime variable (Q n and (t a) ne . 

[00101] Table2B in Figure 3 defines the parameters and variables used in the basic linear 
equations y n versus x n required for preparing the specialized Cartesian graphs in terms of 
pressure and time in column 212; adjusted pseudopressure variable and time in column 213; 
and adjusted pseudopressure and adjusted pseudotime variables in column 214. 

[00102] In rows 204, 205 and 206, the formulas corresponding to coefficients d y (c/), and 
(c 2 ) are given in the case of pressure and time variables in column 212; coefficients (d a \ {c a} \ 
{c a 2\ in the case of adjusted pseudopressure variable and time in column 213; and (d ap ), (c ap i), 
(c ap 2) in the case of adjusted pseudopressure and adjusted pseudotime variables in column 214. 
[00103] In rows 207 and 208, the formulas of the slopes m^and m^M for dual porosity 
reservoir and intercepts b M are given in the case of pressure and time variables in column 212; 
adjusted pseudopressure variable and time in column 213; and adjusted pseudopressure and 
adjusted pseudotime variables in column 214. 

[00104] Figure 4 illustrates three specialized Cartesian graphs of the basic linear 
equations y n versus x n as shown in Tables 2 A and 2B. According to a first series of 
experiment using the same fracture-injection/falloff test data set, the three graphs are three 
straight lines, each having its own slope and intercept. 

[00105] The first series of experiment consists of 21.3 bbl of 2% KC1 water injected at 
5.6 bbl/min over a 3.8 min injection period. In this example, the injection fluid is considered as 
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being a slightly compressible fluid. On the contrary, the reservoir contains a compressible fluid 
that is a dry gas with a gas gravity of 0.63 without significant contaminants at 160°F. 
[00106] The pressure is measured at the surface or near the test interval. The bottomhole 
pressure is calculated from the pressure measurements by correcting the pressure for the depth 
and hydrostatic head. The time interval for each pressure measurement depends on the 
anticipated time to closure. If the induced fracture is expected to close rapidly, pressure is 
recorded at least every second during the shut-in period. If the induced fracture required several 
hours to close, pressure may be recorded every few minutes. The resolution of the pressure 
gauge is very important. The special plotting fiinctions require calculating pressure differences, 
so it is important that a gauge correctly measure the difference from one pressure to the next, 
but the accuracy of each pressure is not critical. For example, consider pressures of 500.00 psi 
and 500.02 psi. The pressure difference is 0.02 psi, so the gauge needs to have resolution on the 
order of 0.01 psi. On the other hand, it doesn't matter if the gauge accuracy is poor. For 
example, if the gauge measures 505.00 and 505.02, then the measurement is within 1% of the 
actual value. Although there is measurement error in the magnitude of the pressure, the 
pressure difference is correct. The analysis is affected by resolution (the difference between 
two measurements), but not necessarily the accuracy. 

[00107] Reservoir pressure is estimated to be approximately 1,800 psi, and the 
bottomhole instantaneous shut-in pressure was 2,928 psi with fracture closure stress observed 
at 2,469 psi. The specialized Cartesian graphs of Figure 4 use the three forms of plotting 
functions defined in the three columns of Tables 2 A and 2B. The method as used in the prior 
art which involves the pressure and time variables evaluates the permeability to be 0.0010 md. 
However, according to the present invention, by using adjusted pseudopressure variable and 
time, the permeability is estimated to be 0.0018 md. And by using adjusted pseudopressure and 
adjusted pseudotime variables, the permeability is estimated to be 0.0023 md. Figure 4 
demonstrates that the fracture- injection/falloff test interpretation is influenced by the pressure- 
dependent properties of the reservoir fluid. Assuming the 0.0023 md permeability estimate is 
correct, then ignoring the pressure-dependent fluid properties by using a pressure and time 
formulation results in a 57% permeability estimate error. 

[00108] According to a second series of experiment, Figure 5 shows three other 
specialized Cartesian graphs of the basic linear equations y n versus x n as defined in Tables 2A 
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and 2B. These three graphs are also represented by three straight lines with different slopes and 
intercepts. 



at 3.3 bbl/min over a 5.2 min injection period. The reservoir contains dry gas with a gas 
gravity of 0.63 without significant contaminants at 160°F. As in the first series of experiment, 
the injection and reservoir fluids are respectively considered as slightly compressible fluid and 
compressible fluid. Reservoir pressure is estimated to be approximately 2,380 psi, and the 
bottomhole instantaneous shut-in pressure was 3,147 psi with fracture closure stress observed 
at 2,783 psi. 

[00110] The specialized Cartesian graph of Figure 5 uses the three forms of plotting 
functions defined in the three columns of Tables 2 A and 2B. The method as used in the prior 
art which involves the pressure and time variables estimates the permeability to be 0.013 md. 
However, according to the present invention, by using adjusted pseudopressure data and time 
as variables, the permeability is estimated to be 0.018 md. By using adjusted pseudopressure 
and adjusted pseudotime variables, the permeability is estimated to be 0.019 md. Once again, 
Figure 5 demonstrates that the fracture-injection/falloff test interpretation is influenced by the 
pressure-dependent properties of the reservoir fluid. Assuming the 0.019 md permeability 
estimate is correct, then ignoring the pressure-dependent fluid properties by using a pressure 
and time formulation results in a 32% permeability estimate error. 

[001 1 1] Both series of experiments also confirm that as pressure approaches and exceeds 
3,000 psi, gas pressure-dependent fluid properties generally will not effect the interpretation 
significantly. However, adjusted pseudovariables are applicable at all pressures and are 
recommended for analyzing all fracture-injection/falloff tests with compressible fluids. 
[00112] Figure 6 illustrates a general flow chart representing a method of iterating the 
measurements and plotting the Cartesian graphs thereof. This graph may apply to the case of 
where the variables are adjusted pseudopressure and adjusted pseudotime. 
[00113] The time at the end of pumping, t ne9 becomes the reference time zero, at step 
600, and the wellbore pressure is measured at At = 0. At steps 602 and 604, calculate the 



[00109] 



The second series of experiment consists of 17.7 bbl of 2% KC1 water injected 



coefficients c a \ 
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[001 14] At step 606, initialize an internal counter n to ne + 1, and test at step 610, if n is 
still below the n max which corresponds to the data point recorded at fracture closure or the last 
recorded data point before induced fracture closure. As is previously said, the time interval for 
each pressure measurement depends on the anticipated time to closure. If the induced fracture 
is expected to close rapidly, pressure is recorded at least every second during the shut-in period. 
If the induced fracture required several hours to close, pressure may be recorded every few 
minutes. 

[00115] If n is below n max , calculate the shut-in time relative to the end of pumping as 
At = t-t ne at step 612. 

[00116] Since the reservoir contains a compressible fluid, its properties will involve the 
calculation of adjusted pseudovariables. At step 614, the adjusted pseudotime variable is 
determined by: 

f (A/) " dAt 

(ta)n = G t V r /)o 7 r- • 111 an embodiment, (tj n is calculated though it is possible to use 

Jo 0*g c t)w 

time as a variable. At step 616, the adjusted pseudopressure variable is determined by: 



P 



(Pw) » pdp 



[00117] At step 622, in Figure 6B, based on the compressibility properties of the 
reservoir fluid, calculate the adjusted pseudopressure variable difference as: 

(Ap a ) n = (Paw)n ~~Par > which can be written as: 

[Pa(P)]„_ x -[Pa(P)] n 



( c t)n 



( t a)n -( t a)n-\ 



[00118] At step 624, calculate the dimensionless before-closure pressure-transient 
adjusted variable (y ap )n defined as: 
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ar 



(Pa)n ~ P 



[00119] At step 626, calculate the dimensionless before-closure pressure-transient 
adjusted variable (x ap ) n defined as: 



( x ap)n - 



c ap\ 



(dap)ne+2 



( d ap)n 



l'/2 



+ £ lVap)j -(dap)j-\]((<a)n -('a)j-\ 
j=ne+3 



( d ap)n 



*nhie j 



c ap20a^n 2 



(d ) W 2 
\ u ap)n l n l ne 



1- 



(*a)/7g+l ^ 
ifa)n j 



1/2' 



[00120] At step 628, increment the internal counter n by 1 and loop back to step 610 to 
test if n is still below n max . 

[00121] At step 610, if n is above m max , Figure 6C indicates that at step 632, prepare a 
graph of (y ap )n versus (x ap ) n . 

[00122] From the graph obtained, and more specifically from the straight line, derive the 
value of the intercept b M which will lead to the evaluation of the reference fracture-face 
resistance R 0 at step 634 using the formula: 



_ 5.615 
R °~ 141.2^(24) r P S f^ 



[00123] However, in order to evaluate the value of the reservoir permeability k, a test at 
step 636 is done in order to determine if the analysis is performed in a dual-porosity reservoir 
system. If it is the case of a single porosity, the value of the slope m M will lead directly to the 
evaluation of the permeability k, at step 640 by calculating the formula as follows, at step 638: 



(141.2)(2)(0.02878)(24) 1 



5.615 



r pSfm M 
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[00124] If it is the case of a dual porosity, the value of a product cok can be evaluated at 
step 650 by calculating the formula as follows, at step 639: 



2 



cok = 



(141.2)(2)(0.02878)(24) 
5.615 



r p S f m M 



[00125] 



Figure 7 illustrates schematically an example of an apparatus located in a drilled 



wellbore to perform the methods of the present invention. Coiled tubing 710 is suspended 
within a casing string 730 with a plurality of isolation packers 740 arranged spaced apart 
around the coiled tubing so that the isolation packers can isolate a target formation 750 and 
provide a seal between the coiled tubing 710 and the casing string 730. These isolation packers 
can be moved downward or upward in order to test the different layers within the wellbore. 

[00126] A suitable hydraulic pump 720 is attached to the coiled tubing in order to inject 
the injection fluid in a reservoir to test for an existing fracture or a new fracture 760. 
Instrumentation for measuring pressure of the reservoir and injected fluids (not shown) or 
transducers are provided. The pump which can be a positive displacement pump is used to 
inject small or large volumes of compressible or slightly compressible fluids containing 
desirable additives for compatibility with the formation at an injection pressure exceeding the 
formation fracture pressure. 

[00127] The data obtained by the measuring instruments are conveniently stored for later 
manipulation and transformation within a computer 726 located on the surface. Those skilled in 
the art will appreciate that the data are transmitted to the surface by any conventional telemetry 
system for storage, manipulation and transformation in the computer 726. The transformed data 
representative of the before and after closure periods of wellbore storage are then plotted and 
viewed on a printer or a screen to detect the slope and the intercept of the graph which may be a 
straight line. The detection of a slope and an intercept enable to evaluate the physical 
parameters of the reservoir and mainly its permeability and face-fracture resistance. 

[00128] The invention, therefore, is well adapted to carry out the objects and to attain the 
ends and advantages mentioned, as well as others inherent therein. While the invention has 



31 



been depicted, described and is defined by reference to exemplary embodiments of the 
invention, such references do not imply a limitation on the invention, and no such limitation is 
to be inferred. The invention is capable of considerable modification, alternation and 
equivalents in form and function, as will occur to those ordinarily skilled in the pertinent arts 
and having the benefit of this disclosure. The depicted and described embodiments of the 
invention are exemplary only, and are not exhaustive of the scope of the invention. 
Consequently, the invention is intended to be limited only by the spirit and scope of the 
appended claims, giving full cognizance to equivalents in all respects. 
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Glossary: 






Ar 




one win$? one face fracture area L 2 ft 2 


br 
°fs 




fracture-face damaee-zone thickness L ft 


°M 




before-closure snpcia1i7pd nlot intprrpnt Himpnsionlpss 

Uvlvlv vlU J HI OUvviailZjvU Ly L lllll/l V/wU I, Ullll^/llOlv/lllV'DO 


B 




formation volume factor dimensionless hhl/STR 


S 


— 


gas formation volume factor, dimensionless, bbl/Mscf 


B 




average £?as formation volume factor dimensionless hhl/lVtscf 


c \ 




V\ O /"* 1 API lt*0 *>t*^P PI IT* A 4"t*0'M PI ^»"% 4" r» 1 -i rni o * /O*** r»V*1 /T 4-3/2 ' 1 /2 1 /2 

oeiore-ciosure pressure-transient analysis vanaoie, twLX , psi cp 


c 2 




lAO'P/XV/a 1 /"X P 1 1 tM*orcni*a ffi«nifl«f nnn1i;f<ip i rn«M ol\1a *"V\^/T 2^.7/2 ♦^o<^^2 -^__l/2 

Deiore-ciosure pressure-transient analysis vanaoie, m /l> t , psi cp 


c a\ 




before-closure pressure-transient analysis adjusted variable, m/Lt 3/2 , psi 1/2 -cp l/2 


c a2 




before-closure pressure-transient analysis adjusted variable, m /L t , psi -cp 


c t 


= 


• • j i 
total compressibility, Lt /m, psi 


C t 


= 


average total compressibility, Lt 2 /m, psi' 1 


d 




before-closure pressure-transient analysis variable, m/Lt 3 , psi/hr 


da 




before-closure pressure-transient analysis adjusted variable, m/Lt , psi/hr 


E 




plane-strain modulus, m/Lt 2 , psi 


u 
n 




JLUIIIlallUIl UllClvIlCbo, 1j, 11 


hf 

J 





fracture height, L, ft 


P 





fracture permeable thickness, L, ft 


J 




index dimensionless 


K 




pcxxiicduiiiiy, , mo 


JU 





dual-norositv bulk- fracture oermeabilitv L 2 md 


T * 




nyaraunc iracture nan lengtn, l, 9 it 


m M 




before-closure specialized plot slope, dimensionless 


n 




index, dimensionless 


P 




pressure, m/i^t , psi 


P 




avpracrp TYrpQciirp m/T t rwi 


Pa 




aUJUalCU piCDoUIC VallaUlC, ill/ LA , pol 


Par 




JiHiiiQfpH tpqptvoi r "nrpQQii rf* vjiti 5iV*1p m/T t^ci 
aujuoi^Li lt-ovi v uii picjouic vcuiduic, iiu i-^i , Uoi 


Paw 


_ 


wellbore adjusted pressure variable, m/Lt 2 , psi 


PaLfD 


= 


dimensionless adjusted pseudopressure variable in a hydraulically fractured 






wpII 


V\A> 


— 


wellbore pressure, m/Lt 2 , psi 


PL f D 




dimensionless pressure in a hydraulically fractured well, dimensionless 


Ap 




pressure difference, m/Lt 2 , psi 






adjusted pressure variable difference, m/Lt 2 , psi 






adjusted pressure variable difference across reservoir zone, m/Lt 2 , psi 
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a) face 


= 


fracture-face adjusted pressure variable difference, m/Lt 2 , psi 


bPcake 


= 


pressure difference across filtercake, m/Lt 2 , psi 


tyface 


= 


pressure difference across fracture-face, m/Lt 2 , psi 




= 


pressure difference across filtrate invaded zone, m/Lt , psi 


*Ppiz 


= 


pressure difference across polymer invaded zone, m/Lt 2 , psi 


&Pres 


= 


pressure difference across reservoir zone, m/Lt 2 , psi 




= 


one wing hydraulic fracture leakoff rate, L 3 /t, bbl/D 


(<lt)g 




one wing hydraulic fracture gas leakoff rate, L /t, bbl/D 


U f 


= 


hydrauhcally fractured well flow rate, L 7t, STB/D 


(U f )g 


= 


hydraulically fractured well flow rate, L 3 /t, STB/D 


r f 




hydraulic fracture radius, L, ft 


r p 

r 


= 


ratio of permeable to gross fracture area, dimensionless 


/ X) 





reference fracture- face resistance m/L 2 t cn/ft 

A A. VA VllV W 11 VAV V VAA W 1UV W A V v A tvil 1 VV Vy ^ AAA/ A -J V ^ VL// A. V 


Re. 





reference fracture- face resistance L 1 ft 1 

A VI WA V1AVW i-1 tiVVVil V 1 VUiUiUllvVk 1 J * A V 


R/s 


= 


fracture-face resistance, L~ l , ft/md 


Rn 





dimensionless fracture-face resistance, dimensionless 


s 


— 


skin, dimensionless 


s f 


= 


fracture-face skin, dimensionless 


s f 


= 


fracture stiffness, m/L 2 t 2 , psi/ft 


t 


= 


time, t, hr 


taLfD 


= 


hydraulically fractured well dimensionless adjusted time, dimensionless 


*n 


= 


time at timestep «, t, hr 






time at the end of an injection, t, hr 


J 


= 


hydraulically fractured well dimensionless time, dimensionless 


' Lne 





fluid volume lost from one wins? of a hvdraulic fracture during an injection T ^ ft"* 


L 




fracture lost width L ft 


x n 




oeiore-ciosure pressure-transient analysis variable, cimensioniess 






i r* i a a. i • j • a i *iii* • 1 

before-closure pressure-transient analysis adjusted variable, dimensionless 


(y a )n 




before-closure pressure-transient analysis adjusted variable, dimensionless 


y n 


= 


before-closure pressure-transient analysis variable, dimensionless 


z 


— 


gas deviation factor, dimensionless 


z 




average gas deviation factor, dimensionless 


Greek 






Li 




viscosity, m/Lt, cp 


M 




average viscosity, m/Lt, cp 






gas viscosity, m/Lt, cp 


* 




porosity, dimensionless 


0) 




natural fracture storativity ratio, dimensionless 



34 



